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Abundance Patterns in the Interstellar Medium of the SO Galaxy 
NGC 1316 (Fornax A) Revealed with Suzaku 
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Abstract 

The Suzaku X-ray satellite observed the nearby SO galaxy NGC 1316, a merger remnant aged 3 Gyr. 
The total good exposure time was 48.7 ksec. The spectra were well represented by a two-temperature 
thermal model for the interstellar medium (ISM) plus a power-law model. The cool and hot temperatures 
of the thermal model were 0.48±0.03 and 0.92±0.04 keV, respectively. The excellent spectral sensitivity 
of Suzaku enables for the first time to measure the metal abundances of O, Ne, Mg, Si, and Fe in the ISM. 
The resultant abundance pattern of O, Ne, Mg, Si, and Fe is close to that of the new solar abundance 
determined by Lodders (2003). The measured abundance pattern is compared with those of elliptical 
galaxies and an SO galaxy, observed with Suzaku. Considering the metal-enrichment from present Type la 
supernovae, the near-solar abundance pattern of the ISM in NGC 1316 indicates an enhanced a/Fe ratio 
of stellar materials in the entire galaxy, like in giant elliptical galaxies. 

Key words: galaxies: individual (NGC 1316), galaxies: interstellar medium, galaxies: abundances 
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1. Introduction 

Hot X-ray emitting interstellar medium (ISM) carries 
important information about the history of star forma- 
tion and evolution of galaxies. The metals in the ISM in 
early-type galaxies have been enriched by type la super- 
novae (SNe la) and stellar mass loss. Since SNe la do not 
produce O or Mg, abundances of these elements in the 
ISM reflect stellar abundances. Thus, the abundance pat- 
tern of the ISM can play a key role in studying the origin 
of metals. 

The ASCA satellite first enabled us to measure the 
metal abundances in the ISM of the early-type galax- 
ies (e.g., Awaki et al. 1994; Loewenstein et al. 1994; 
Mushotzky et al. 1994; Matsushita et al. 1994). These 
galaxies have about 1 solar iron abundance, employ- 
ing a multi-temperature plasma model (Buotc & Fabian 
1998), or considering uncertainties in the Fe-L atomic data 
(Arimoto et al. 1997; Matsushita et al. 1997; Matsushita 
et al. 2000). Chandra and XMM- Newton measured the 
metal abundances and spatial distributions of metals in 
the ISM of the early-type galaxies, because of their large 
effective area and good angular resolution (Humphrey et 
al. 2004; Humphrey & Buote 2006; Werner et al. 2006; 
Tozuka & Fukazawa 2008). In particular NGC 5044, 
NGC 1399, and NGC 4636, the abundances of non- 
Fe elements have good constraints with Chandra and 
XMM (Buote et al. 2003; Ji et al. 2009). The observed 
abundance patterns are consistent with the solar ratio. 



However, O and Mg abundances of relative faint X-ray 
galaxies still remain highly uncertain, due to a highly 
asymmetric energy response in the low-energy region be- 
low 1 keV and a strong instrumental Al line of these 
missions, respectively The Suzaku XIS (Koyama et al. 
2007) has an improved line spread function due to a very 
small low-pulse-height tail below 1 keV, with a lower and 
more stable background level compared to the Chandra 
ACIS and the XMM-Newton EPIC. The abundance pat- 
terns in the ISM of several giant elliptical galaxies, includ- 
ing NGC 1399 and NGC 1404 (Matsushita et al. 2007), 
NGC 720 (Tawara et al. 2008), and NGC 4636 (Hayashi 
et al. 2009), have been revealed by recent Suzaku observa- 
tions. The derived Fe abundances of the ISM are about 1 
solar, and the abundance patterns are close to that of the 
new solar abundance determined by Lodders (2003). From 
these results, the contribution of SNe la and stellar mass 
loss to the metal enrichment of the ISM and intracluster 
medium were discussed. However, NGC 4382 is the only 
SO galaxy whose ISM abundance pattern has so far been 
derived with Suzaku (Nagino & Matsushita 2010). The 
O /Fe ratio of the ISM in this galaxy is smaller by a fac- 
tor of two as compared to than those in the four elliptical 
galaxies. 

In terms of the cosmology dominated by cold dark mat- 
ter, it is considered that at least some elliptical galax- 
ies were formed by a hierarchical formation scenario, in 
which larger spheroidals were assembled relatively late 
through mergers of the late-type galaxies of comparable 
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mass (De Lucia et al. 2006 and references therein). In con- 
trast, SO galaxies in clusters of galaxies are considered to 
be formed through morphological transformations, from 
infalling gas-rich spiral galaxies to gas-poor SO galaxies. 
This is because spiral galaxies were much more common, 
and SO galaxies were much rarer, in distant than nearby 
galaxy clusters (e.g. Poggianti et al. 2009; Desai et al. 
2007; Smith et al. 2005; Postman et al. 2005). In spite of 
these differences, the elliptical and SO galaxies have been 
treated together in most of the past studies. However, 
since there are possible differences in the formation sce- 
nario, we need to study the elliptical and SO galaxies sep- 
arately. 

The abundance pattern of stars reflects their forma- 
tion history, since a longer formation time provides a 
higher concentration of trapped SNe la products to the 
stars. If the elliptical and SO galaxies have differ- 
ent formation histories, their stellar metallicity and/or 
present SN la rate may be different. The stellar metal- 
licity of the early-type galaxies has often been investi- 
gated in optical observations of their central regions (e.g., 
Thomas et al. 2005; Bedregal et al. 2008; Walcher et 
al. 2009). However, there may be systematic uncertain- 
ties in the assumption of the age-population of stars, 
and in atomic physics. Also the values of solar abun- 
dances have changed enormously, considering the three- 
dimensional hydrostatic atmosphere and non-local ther- 
modynamic equilibrium (Asplund 2005). In contrast, X- 
ray observations of the ISM, show much simpler atomic 
data and temperature structures as compared to those in 
the optical spectra, and we can observe the entire region 
of each galaxy, although the metals in the ISM are a mix- 
ture of those from stars and recent SNe la. Considering 
these, further X-ray studies of SO galaxies are thought to 
be important. 

NGC 1316 is a peculiar SO galaxy, with numerous tidal 
tails. This galaxy has a surprisingly high central surface 
brightness, and a lower velocity dispersion (a = 227±33 
km s -1 ; Goudfrooij et al. 2001) for its galaxy luminos- 
ity than other elliptical galaxies with similar luminosi- 
ties (e.g., Schweizer 1980; Schweizer 1981; D'Onofrio et 
al. 1997). Due to these features, NGC 1316 has been 
thought to have undergone a major merger relatively re- 
cently. The merger age is estimated to be ^3 Gyr, based 
on the age of the bright globular clusters (Goudfrooij et 
al. 2001). Therefore, NGC 1316 is good target for inves- 
tigation of the galaxy formation process through mergers. 
Measurements of the absorption-line indices of the central 
region, or 1/10 r e of NGC 1316, indicates a small over- 
abundance of a-to-Fe ratio of [a/Fe]=0.15 (Thomas et 
al. 2005), which is smaller than the [a/Fe]~0.3, in center 
of the giant elliptical galaxies (e.g. Thomas et al. 2005; 
Pipino et al. 2009). Here, r e is the effective radius of 
the galaxy. These line indices also indicate a somewhat 
younger age of 3-5 Gyr (Bedregal et al. 2008). 

ROSAT observations of NGC 1316 detected a thermal 
emission component from hot ISM (Kim et al. 1998), using 
Chandra data of an good exposure of 25 ksec, Kim & 
Fabbiano (2003) found that the ISM has a temperature 
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Fig. 1. The 0.4-5.0 kcV XIS image of NGC 1316, overlaid on 
an optical contour map taken by Digitized Sky Survey. The 
observed XIS0, 1, and 3 images were added on the sky co- 
ordinate, and smoothed with a a = 5.2 arcsec Gaussian pro- 
file. Estimated contributions of the cosmic X-ray background 
(CXB) and the instrumental non-X-ray background (NXB) 
were subtracted, although vignetting was not corrected. The 
region used for spectrum extraction is indicated by a green cir- 
cle. NGC 1317, expressed by a reg circle, was excluded from 
the background and galaxy spectra. 

in the 0.5-0.6 keV range and a 0.3-8 keV luminosity of 
3.1 xlO 40 erg s _1 . The metal abundance of the ISM were 
reported to be 0.2-1.3 solar adopting the solar abundance 
table by Anders & Grevcssc (1989). NGC 1316 hosts radio 
lobes, that generate inverse- Compton X-rays as well as 
synchrotron radio emission (e.g. Feigelson et al. 1995; 
Kaneda et al. 1995; Isobe et al. 2006). Although this 
means that NGC 1316 harbors an active galactic nucleus 
(AGN), Iyomoto et al. (1998) and Kim & Fabbiano (2003) 
showed that the present activity of the AGN is very low. A 
recent Suzaku observation detected the inverse- Compton 
X-ray up to 20 kcV utilizing Suzaku from the west lobe 
(Tashiro et al. 2009). 

In this paper, we analyze Suzaku XIS data of 
NGC 1316. We use a distance of 18.6 Mpc to NGC 1316 
based on the Hubble Space Telescope measurements of 
Cepheid variables (Madore et al. 1999), and a redshift of 
0.005871 including proper motion (Longhetti et al. 1998). 
We use the new solar abundances in Lodders (2003). 
Unless noted otherwise, the quoted errors are refer to 90% 
confidence interval for a single interesting parameter. 

2. Observation and Data Reduction 

Suzaku observed NGC 1316 in 2006 December, us- 
ing the XIS (Koyama et al. 2007). The XIS consists of 
two front-illuminated (FI: XIS0 and XIS3) CCD cam- 
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Table 1. Best-fit parameters for the apec components + power-law models.* 
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Parameters 



(i) 



(ii) 



OTHERS 



iVormoTHERS ' 

kTiGM 

Abundance 
NormicM * 

^Tmwh 
Abundance 

NormMwn t 

kT L HB 

Abundance 
Norm L nB t 

X 2 /d.o.f. 



(cm 2 ) 

(keV) 
(solar) 

(keV) 
(solar) 

(keV) 
(solar) 



1.6 (fix) 
1.10±0.07 
1.92 xlO 20 (fix) 
n 74+0-06 

u - ' 4 -0.05 

0.3 (fix) 
1.88±0.36 

30+ - 03 

u - ou -0.04 

1 (fix) 
0.77±0.17 

0.1 (fix) 
1 (fix) 



2.79 



+0.59 
-0.58 



490/359 



1.6 (fix) 
1.10±0.07 
1.92 xlO 20 (fix) 
n 74+0 .06 

u -'^-0.05 

0.3 (fix) 
1.91±0.35 

30+ 03 

1 (fix) 
0.68±0.15 

0.1 (fix) 
1 (fix) 
2.73±0.60 

490/359 



* The apec components for spectra in the background region of NGC 1316 with the ICM component of the Fornax cluster, absorbed or 
non-absorbed MWH component, LHB component for the Galactic emission, and a power-law model for CXB. 

* Normalization of the power-law component divided by the solid angle, Q u , assumed in the uniform-sky ARF calculation (20' radius), in 
units of 10 -3 0° photons keV -1 cm -2 s _1 arcmin -2 at 1 keV. 

i Normalization of the apec components divided by the solid angle same as the normalization of apec, Norm = J n c nndV / (47r (1 + 
z) 2 D^J /Q" XlO - 17 cm -5 arcmin -2 , where D\ is the angular distance to the source. 



eras and one back-illuminated (BI: XIS1) CCD camera. 
The averaged pointing direction of the XIS was at (a, 
5)=(3 h 22 m 40?4, -37°12'10."4), right on the NGC 1316 
nucleus. The XIS was operated in normal clocking mode 
(8 s exposure per frame), with the standard 5x5 and 
3x3 editing mode. We processed the XIS data using 
the "xispi" and "makepi" ftool tasks and CALDB files of 
2008-08-25 version. Then, the XIS data were cleaned by 
assuming thresholds on the Earth elevation angle of > 5° 
and the Day-Earth elevation angle of > 20°. We also dis- 
carded data taken when the time after a spacecraft exit 
from the south Atlantic anomaly is less than 436 sec. We 
created a 0.5-2 keV light curve for each sensor, with 540 
sec binning, to exclude periods of anomalous event rates 
higher or less than ±3a from the mean. After this screen- 
ing, the remaining good exposures were 48.7 ksec for both 
FIs and BI. Event screening with cut-off rigidity was not 
performed. 

The spectral analysis was performed with HEAsoft 
version 6.5 and XSPEC 12.4. To subtract the non- 
X-ray background (NXB), we employed the Dark-Earth 
database using the "xisnxbgen" ftool task (Tawa et al. 
2008). We generated two different ancillary response files 
(ARFs) for the spectrum of each region; one assumed uni- 
form sky emission, while the other utilized the observed 
XIS1 image by the "xissimarfgen" ftool task (Ishisaki et 
al. 2007). In the ARFs, we also included the effect of 
contamination on the optical blocking filters of the XIS. 



3. Analysis and Results 

3.1. Extraction Regions of Spectra 

The observed X-ray image in the 0.4-5 keV range is 
shown in figure 1, superposed on an optical contours from 
the Digitized Sky Survey. We used the following two re- 
gions to extract spectra; one is NGC 1316 region, which is 
a circular region with a 5' radius shown in a green circle 
in figure 1. The other is background region, which is the 
entire XIS field of view excluding the NGC 1316 region. 
We excluded the region around NGC 1317, shown in a red 
circle in figure 1. In addition, three calibration sources, 
with emission peaks at 5.9 keV, are located at the corners 
of the XIS. We included these regions to improve the pho- 
ton statistics because the ISM emits almost no photons 
above 5 keV. 

3.2. Estimation of Background Spectra 

In order to estimate the background components, we 
first fitted the spectra extracted over the background re- 
gion. Since the NXB component was already subtracted, 
the remaining background consisted of extra-galactic cos- 
mic X-ray background (CXB), and emission from our 
Galaxy. We assumed a power-law model for the CXB com- 
ponent, and a two temperature model from the Galactic 
emission. Empirically, one thermal component represents 
the sum of solar wind charge exchange (SWCX) and local 
hot bubble (LHB), while Milky Way halo (MWH) emits 
the other thermal component (Yoshino et al. 2009). We 
also added a thin thermal plasma emission from intra- 
cluster medium (ICM) of the Fornax cluster (Tashiro et 
al. 2009). We thus, fitted the spectra with the follow- 
ing model : phabs x (power-law + apecicM + apecMWH) 



S. Konami et al. [Vol. 




Energy (keV) Energy (keV) 

Fig. 2. NXB-subtracted XISO (black) and XIS1 (red) spectra of NGC 1316 (left column), and those of the background region 
(right column), shown without removing the instrumental responses. The top two and bottom two panels employ the one- and 
two-temperature model for the ISM, respectively. Black and red lines show the best-fit model for the XISO and XIS1, respectively. 
For simplicity, only the model components for XIS1 spectra are shown. Orange and purple lines are the ISM components, cyan is 
emission from the ICM of Fornax cluster by apecicM > magenta and yellow are the Galactic background emission by apecMWH and 
apecLHB i while blue and green are the CXB and (LMXB + radio lobe) components, respectively. The background components, except 
power-lawoTHERSi are common between the on-source and background spectra, but scaled to the respective data accumulation area. 
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Table 2. Summary of the best- fit parameters for the NGC 13f6 and the background region with IT for ISM model, 2T for ISM 
model, and 2T for ISM adding 10% systematic error of Fe-L model. 



Parameters 




IT for ISM 


2T for ISM 


Fe-L 


kTi T 


(kcV) 


0.62±0.01 


0.48±0.03 


40 +IU)8 

"■^"-0.03 


Normn * 




2 10+ 017 

z ' lu -0.18 


1 04+ 21 


95+ 16 


kT 2 T 


(kcV) 




0.92±0.04 


OQ+0 05 
U - JO -0.04 


Norm2T * 




- 


0.80jJiS 


70+ 12 

u - ' u -0.22 





(solar) 


n 7i+ - 41 

u - ' i -0.23 


0.79i ji 


n 71+0.50 

u - ' ± -0.26 


Ne 


(solar) 


1 O 1 T^U.O < 

i -° 1 -0.35 


0.63t°;g 


0-66±g;I| 


Ms, Al 


(solar) 


0.43 + n?? 
^ ^ — 0.13 


w ^^ — 0.23 


' —0.33 


Si S Ar Ca 


[ solar"! 


31+° ?° 

U - J -0.12 


48 + S o= 

"■^°-0.38 


64 + nne 

"■ u -0.25 


Fp Ni 


1 Qnl £1 T 1 

i auicii j 


u -^-0.09 


u - '°-0.15 


1 01+ 075 


O/Fe 




1 61+° I 3 , 


1 08 + °1? 


70+°^ 

u " ' u -0.25 


Ne/Fe 


( Qnl ar 1 

l oUlal J 


q nn-|-n 57 


86+ 68 


65+ 083 

u - ud -0.52 


Mg/Fe 


(solar) 


0.98±°if 


0.75t° 2 25 


7fi+ - 26 
u - ' u -0 06 


Si/Fe 


(solar) 


f) 70 +a3 ° 

U. /U_ 27 


fifi +0 - 21 

U.DD_ 23 


U.DO_ 20 




(keV) 


n 77+O.O8 
u - ' ' -0.06 


n 79+O.O5 
u -' z -0.04 




AormicM * 




2.51±0.51 


2.83±0.36 






(keV) 


0.30±0.03 


n 90+0.02 

U - zS -0.03 


0.28±0.02 


Aorm M \VH * 




1 44+ - 27 


1 2fi +0 - 23 

l.^D_ 2 8 


1 O o+0.21 

l- z o_ 14 




(keV) 


0.1 (fix) 


0.1 (fix) 


0.1 (fix) 


Aor7TiLHB * 




4.52±0.98 


4 50+ 1 - 00 

*- ou -1.40 


4 Q7+1-41 


X 2 /d.o.f. 




970/703 


921/701 


859/692 



* Normalization of the vapec and apec components 
S OURCE_RATIO.REG J / [An {l + z) 2 D 2] xlQ -17 



scaled with a factor of S0UR.CE_RATI0_REG / AREA, which 
cm -5 arcmin~ 2 , where Da is the angular distance to the source. 



is Norm 



+ apecLHB- Here, the "phabs" factor represents the 
photoelectric-absorption, whose column density was fixed 
to the Galactic value of 1.92 x 10 20 cm~ 2 in the direction of 
NGC 1316. The "apecicivi", "apecMWH" and "apecLHB" 
terms mean the thermal emission, in terms of APEC 
plasma model (Smith et al. 2001), from the ICM, MWH, 
and LHB (and SWCX), respectively. We assumed that the 
apccicM has a metal abundance fixed at 0.3 solar, while 
apecMWH and apecLHB 1-0 solar. The redshift of apecicM 
were also fixed at 0.005871, while apecMWH and apecLHB 
zero. The power-law slope of the CXB emission was fixed 
at T = 1.4 (Kushino et al. 2002). The spectra from the 
BI and FI CCDs were simultaneously fitted in the 0.4-5.0 
and 0.5-5.0 keV range, respectively. 

This model well reproduced the spectra, with x 2 /d.o.f. 
= 455/360. However, the strength of the power-law rep- 
resenting the CXB was about twice as high as that of 
Kushino et al. (2002). We consider that this is due to con- 
tributions by other sources in NGC 1316, such as the radio 
lobes and integrated low mass X-ray binaries (LMXBs). 
The index of the power-law component from the radio- 
lobes is 1.68 (Tashiro et al. 2009). The integrated spectra 
from discrete sources in the early type galaxies are ap- 
proximated with a power-law model with an index of 1.6 
(e.g., Xu et al. 2005, Randall et al. 2006). Therefore, 
we added a power-law component, "power-lawoTHERs" j 
with its index fixed at 1.6, and then fitted the spectra 



with the following model (i): phabs x (power-lawcxB + 
power-lawoTHERs + apecicM + apec M wii) + apec LH B- 
The index and normalization of power-lawcxB were fixed 
at the values by Kushino et al. (2002). The results of the 
fit are summarized in table 1. 

We also examined the following model (ii): phabs 
x (power-lawcxB + power-law THERS + apecicM) + 
apecMWH + apecLHB- The difference of this model from 
model (i) is whether or not the MWH component is ab- 
sorbed by the Galactic column (e.g., Yamasaki et al. 2009; 
Konami et al. 2009). We summarize the results in table 1. 
There was no significant difference in the derived param- 
eters between model (i) and (ii). Hereafter, we adopted 
model (i) as the background model for the spectral fitting, 
unless otherwise stated. 

To evaluate the Galactic emission, we also fitted the 
O VII and O VIII lines of the spectra of the background 
region with Gaussians and a power-law model in the 0.5- 
0.7 keV range. The resultant intensities of the OVII and 
O VIII, without correcting for absorption, are 9.2± 5 4 7 3 and 
2.4lgg photons cm -2 s _1 sr _1 , respectively. These values 
are consistent with the previously reported values from 
other sky region (McCammon et al. 2002; Sato et al. 2007; 
Yamasaki et al. 2009; Konami et al. 2009). Thus, the 
Galactic foreground emission in the present sky direction 
is considered typical. 
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Fig. 3. Confidence contours between the metal (O, Ne, Mg, and Si) and Fe abundances, determined with the two-tem- 
perature model for the ISM. The black, red, and green contours represent 68%, 90%, and 99% confidence ranges, respec- 
tively. The red dashed contours indicate 90% confidence range determined with the one-temperature model for the ISM. 



3.3. Emission from NGC 1316 

We fitted the spectra of the NGC 1316 region with the 
model: phabsxvapccigM + background models. In this 
model, "phabs" represents the Galactic absorption in the 
direction of NGC 1316, fixed at A H = 1.92 x 10 20 cnr 2 . 

The ISM emission of NGC 1316, vapecisM, was repre- 
sented by one-temperature (IT) or two-temperature (2T) 
models, employing the vapec code (Smith et al. 2001). 
The abundances of He, C, and N were fixed to the solar 
value. We also divided the other metals into five groups 
as O, Ne, (Mg & Al), (Si, S, Ar, Ca), and (Fe & Ni), 
based on the metal synthesis mechanism of SNe, and al- 
lowed each group to vary. In the 2T modeling, the two 
vapec components were assumed to share the same ele- 
mental abundances. Using this model, we simultaneously 
fitted the spectra of the NGC 1316 and background re- 
gions using the energy ranges of 0.4-5.0 keV for the BI 
CCD and 0.5-5.0 keV for the FI CCDs. The normaliza- 
tion of power-lawoTHERS w &s not linked between the two 



regions, because the number of LMXBs may increase to- 
ward the galactic center, and the brightness of the radio- 
lobes may be different. All the other background compo- 
nents were constrained to be the same between the two 
regions. The normalizations of the ISM components in 
the background region were fixed at 0. 

The results of these fits are summarized in figure 2 
and table 2. Several emission lines seen around 0.5-0.6 
keV, 0.6-0.7 keV, ~ 1.3 keV, and - 1.8 keV are iden- 
tified with Ka lines of OVII, OVIII, MgXI, and Si XIII, 
respectively. The emission bump around 0.7-1 keV corre- 
sponds to Fe-L complex, as well as to K-lines from Ne IX 
and NeX. These fits are not formally acceptable, mainly, 
because of discrepancies between the FI and BI spectra 
of the background region above 2 keV. Although this dis- 
crepancy might be related to some remaining problems 
with the NXB subtraction, they are unlikely to affected 
the ISM modeling. The IT model results in kT = 0.62 
keV. The fit statistics shown in table 2 clearly favor the 
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2T model, which gives the two temperatures as 0.48 and 
0.92 kcV. The abundance of Fe does not depend very 
much on the temperature structure of the ISM. Adoption 
of the 2T model shows, the Fe abundance of the ISM as 
0.58-1.15 solar. This value is consistent with the ISM 
metallicity of 0.4-2.1 solar derived from the Chandra ob- 
servation derived by Kim & Fabbiano (2003). Here, the 
value was converted using the solar abundance table by 
Lodders (2003). 

In order to examine the abundance ratios rather than 
their absolute values, we calculated confidence contours 
between the abundances of metals (O, Ne, Mg, and Si) 
and of Fe, using the 2T model for the ISM. The results 
are shown in figure 3 and table 2. The elongated shape of 
the confidence contours indicates that the relative values 
are determined more accurately than the absolute values. 
Figure 4 the shows metal-to-Fe ratios in NGC 1316 using 
the IT and 2T models for the ISM, which were derived 
from the two-parameter confidence contours shown in fig- 
ure 3. The O/Fe ratio varies from 1.7 with the IT model 
to 1.1 with the 2T model in solar units. On the other 
hand, the Ne/Fe ratio changes considerably, from —3 so- 
lar with the IT model to ~1 solar with the 2T case. We 
note that the Ne abundance is not reliably determined 
due to an overlap with the strong and complex Fe-L line 
emissions. The IT and 2T models give similar values for 
the Mg/Fe and Si/Fe ratios. 

3.4- Systematic Uncertainties in the Abundance Ratios 

We also tried to fit the same spectra with the three- 
temperature vapec model for the ISM. The derived tem- 
perature of the third component has a large error, and 
the Fe abundance increased to 1.03to'29 m solar units. 
Meanwhile, the derived abundance ratios are consistent 
with those from the 2T model. To examine the depen- 
dence on the power-lawoTHERS component, we fitted the 
spectra and allowed the index to vary freely. Then, the in- 
dex decreased to 1.20, but the abundances did not change 
significantly. Systematic uncertainty associated with the 
contaminants on the XIS optical blocking filter is less than 
the statistical errors. 

The derived ICM component in subsection 3.3 are in- 
cluding both ICM and ISM components. Considering the 
point spread function of Suzaku X-ray telescope, ~ 30 % 
of the ICM component in the background region escaped 
from the ISM emission in the galaxy region (Serlemitsos 
et al. 2007). In addition, the ISM emission may extend 
up to the background region. Therefore, we added an- 
other temperature component in the background spectra, 
and fitted the spectra again. We got almost the same val- 
ues for the ISM component in the galaxy region. This is 
because the normalization of the ICM component in the 
galaxy region is only 10 % of that of the ISM, as shown 
in figure 2. 

Because the fits were not formally acceptable, we ex- 
plored using the C-statistic giving less biased values than 
X 2 statistic for Poisson-distributed data (Humphrey ct al. 
2009). The results of C-statistics fitting hardly improved, 
X 2 /d.o.f.=939/701, and the parameters and their errors 



were almost the same. The large x 2 mostly comes from 
the systematic errors of XIS. For examples, because the 
response matrix of XISs around the Si edge has some prob- 
lems, the energy range between 1.82 and 1.84 keV was 
also ignored in the spectral fitting (Koyama et al. 2007). 
Then, the x 2 /d.o.f improved to 906/691. Furthermore, it 
is known the discrepancies between the FI and BI spec- 
tra of the background region above 2 keV. Ignoring this 
energy range in the background spectra, the x 2 /d.o.f. be- 
come 666/578, and the parameters were almost the same 
within a few %. The large x 2 are almost due to systematic 
uncertainties in the background and calibrations of XIS. 

We also investigated the systematic uncertainties in 
the atomic physics of the Fe-L lines as pointed out in 
Matsushita et al. (2000) and Matsushita et al. (2007). 
Convolving the response matrix of XIS, the difference in 
the Fe-L lines between APEC (Smith ct al. 2001) and 
MEKAL (Mewe et al. 1985) models are typically 10%. 
Therefore, we fitted spectra in the same way as in section 
3.3 but including 10 % systematic errors within an energy 
range of 0.75-1.2 keV. The derived parameters are sum- 
marized in table2, and x 2 /d.o.f became 859/692. The er- 
ror range in absolute value of each element increased, and 
the derived Fe abundance became 0.7-1.8 solar. However, 
increases in the error ranges of the abundance ratios were 
smaller. Furthermore, ignoring an energy range of 1.82- 
1.84 keV, and above 2 keV in the background spectra, the 
X 2 /d.o.f became 605/569. 

4. Discussion 

We successfully measured the abundance patterns of 
O, Ne, Mg, Si, and Fe in the ISM of the SO galaxy 
NGC 1316 with Suzaku. Figure 4 summarizes the results. 
We used the new solar abundance table by Lodders (2003). 
Adopting the result of the 2T model fit, the abundance 
pattern in the ISM in NGC 1316 is shown to be close to 
that of the new solar abundance. The abundance patterns 
of the SN II and SN la yields are also plotted in figure 4. 
Here, the SN II yields by Nomoto ct al. (2006) refer to an 
average over the Salpeter initial mass function of stellar 
masses from 10 to 50 Mq, with a progenitor metallicity 
of Z = 0.02. The SNe la yields were taken from the W7 
model Iwamoto et al. (1999). The abundance pattern in 
the ISM of NGC 1316 lies between those of SN II and SN 
la, and is consistent with their mixing. Assuming that 
the O/Fe ratio of the SN II yield is 3.5 solar (Nomoto et 
al. 2006), then -70% of Fe in the ISM is synthesized by 
SNe la, and originates from present SNe la and from those 
trapped in stars. Since the 2T and 3T models give the Fe 
abundance of the ISM of 0.6-1.6 solar, the Fe abundance 
from SNe la would be 0.4-1.1 solar. 

The Fe abundance synthesized by present SNe la in 
an early-type galaxy is proportional to A'/g^6>sN/a* (see 
Matsushita et al. 2003 and Nagino & Matsushita 2010 for 
details). Here, is the Fe mass synthesized in one 

SN la, #sn is the SN la rate, and a* is the stellar mass 
loss rate. Fe in the hot ISM of the galaxy is mainly pro- 
duced from SNe la, since SNe II synthesized far more O 
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Fig. 4. Abundance ratios of O, Ne, Mg, and Si to Fe for the one-temperature (red) and two-temperature (blue) model 
of the ISM in NGC 1316. Abundance patterns of NGC 4382 (purple; Nagino & Matsushita 2010), NGC 1332 (orange; 
Humphrey & Buote 2006), NGC 720 (black crosses; Tawara et al. 2008), NGC 1399 (black closed squares; Matsushita et 
al. 2007), NGC 1404 (black open squares; Matsushita et al. 2007), NGC 4636 (black stars; Hayashi et al. 2009), and 
NGC 5044 (black diamond; Buote et al. 2003 and Komiyama et al. 2009) are also shown. Solid and dot-dashed lines rep- 
resent the number ratios of metals to Fe for the SN II and SN la products (Iwamoto et al. 1999; Nomoto et al. 2006). 



and Mg than Fe. We used the mass-loss rate from Ciotti 
et al. (1991) and assumed the age to be 13 Gyr, which 
is approximated by 1.5 x 10 _11 LBi 1 5 L ' 3 M Q yr _1 , where ti$ 
is the age in units of 15 Gyr and Lb is the B-band lu- 
minosity. Mp e produced by one SN la explosion is likely 
to be ~O.6M (Iwamoto et al. 1999). Adopting 0.1-0.5 
SN Ia/100yr/ 10 10 i B as the optically observed SN la rate 
(e.g., Mannucci et al. 2008 and references therein , Blanc 
et al. 2004, Hardin et al. 2000, and Cappellaro et al. 1997), 
the resultant Fe abundance is 2.9-14.5 solar, when con- 
sidering only the SNe la contribution. Therefore, the ex- 
pected SNe la contribution is sufficiently high to supply 
the Fe in the ISM of this galaxy. 

As indicated from the optical line indices (Thomas et 
al. 2005; Bedregal et al. 2006; Bedregal et al. 2008), 
NGC 1316 may contain younger stellar populations. In 
this case, both the stellar mass loss rate and SN la rate 



might be higher (Ciotti et al. 1991; Ciotti & Ostriker 
2007). The evolution of SN la rate is approximated by 
power- law model: oc t~ ls L B (Ciotti & Ostriker 2007; 
Greggio 2005), where t is the age of stars. Then, assum- 
ing an age of 3 Gyr, the contribution of SNe la to the Fe 
abundance decreases by a factor of 0.7. Since a significant 
proportion of stars in NGC 1316 might be older than 3 
Gyr, the effect of stellar age on the metal-enrichment by 
present SNe la would be less than a few tens of a per- 
cent. Then, the difference in the abundance pattern of 
the ISM of the early-type galaxies reflects the difference 
in the abundance pattern of the stars in the entire galaxy. 

Since a longer star-formation time-scale yields more 
SNe la products in stars, the differences in the a/Fe 
ratio in stars can constrain the star-formation histories. 
The abundance pattern of O, Ne, Mg, Si, and Fe in the 
ISM of NGC 1316 is close to that of the elliptical galax- 
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ies, have been obtained by Chandra, XMM, and Suzaku 
(e.g.Humphrey & Buote 2006, Werner et al. 2006, and 
Matsushita et al. 2007). We show the abundance pat- 
terns of NGC 4382, NGC 1332, NGC 720, NGC 1399, 
NGC 1404, NGC 4636, and NGC 5044 researched previ- 
ously (Nagino & Matsushita 2010, Humphrey et al. 2004, 
Tawara et al. 2008, Matsushita et al. 2007, Hayashi et al. 
2009, Buote et al. 2003, and Komiyama et al. 2009) in 
figure 4. Considering the expected high contribution from 
the present SNe la, the observed solar abundance pattern 
in the ISM of NGC 1316 indicates an overabundance of 
[a/Fe] in stars in the entire galaxy as in these elliptical 
galaxies, although optical observations indicate that stars 
in the central 1/10 r e region of NGC 1316 contain more 
SN la products than those in the giant elliptical galaxies 
(Thomas et al. 2005). Suzaku observed the ISM of the 
entire galaxy and revealed the abundance pattern in stars 
over the entire galaxy. 

The abundance pattern of the SO galaxy, NGC 4382 
was also investigated with Suzaku (Nagino & Matsushita 
2010). NGC 4382 is a normal SO galaxy located on the 
outskirts of the Virgo cluster. Figure 4 also shows that 
the O/Fe ratio of the ISM in NGC 4382 is smaller, by 
a factor of two, as compared to that of NGC 1316. The 
different abundance patterns of these two SO galaxies re- 
flect a higher amount of SNe la products in the ISM of 
NGC 4382. Measurements of optical line indices indicate 
that stars in the most central region (< l/8r e ) of these 
two galaxies have similar [a/Fe] ratios (Thomas et al. 
2005; McDermid et al. 2006). As shown above, the ef- 
fect of stellar age on the Fe abundance of the ISM from 
the present SNe la might be small, but stars in the en- 
tire NGC 4382 may contain more SNe la products than 
those in the entire NGC 1316. These results suggest dif- 
ferences in the metal enrichment histories and hence in 
the formation and evolution histories of the ISM in the 
merger remnant, NGC 1316, and the normal SO galaxy, 
NGC 4382. Optical observations reveal that the a/Fe ra- 
tio of the stars in the central regions of elliptical and SO 
galaxies depends on the system mass (e.g., Thomas et al. 
2005; Bedregal et al. 2008). The temperature of the ISM 
of NGC 4382 is about 0.3 kcV (Nagino & Matsushita 2010) 
while and those of NGC 1316 and the ellipticals are more 
than 0.5 keV, as shown in figure 4. Since the tempera- 
ture of the ISM reflects the system mass, the differences 
in the abundance patterns in the ISM of the two galax- 
ies may also reflect the dependence of the a/Fe ratio of 
stars in the entire galaxy on the system mass. In addition, 
the abundance pattern of another SO galaxy, NGC 1332 
was searched with XMM and Chandra (Humphrey et al. 
2004). Although the temperature of the ISM of NGC 1332 
is > 0.5 keV, the a/Fe has large error to compare with 
other SO galaxies, as shown in figure 4. Further studies 
to obtain, more samples of the elliptical and SO galaxies 
with different ISM temperatures would enable us to un- 
derstand the formation process of the early-type galaxies 
in more detail. 
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